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1 .  INTRODUCTION 

Tnis  presentation  constitutes  the  final  report  to  the  Air  Force 
Geophysical  Laboratory,  LYP  under  Contract  No.  ESD  85-609.  In  this  report  we 
present  a  description  of  the  three-dimensional  NOAA/LAMP  mesoscale  prediction 
model  which  contains  explicit  predictions  of  the  evolving  fields  of  cloud  and 
rain  water.  We  show  some  examples  of  numerical  simulations  made  using  the 
model,  and  discuss  the  success  and  limitations  of  the  present  version.  We 
then  outline  the  status  of  other  efforts  designed  to  overcome  deficiencies 
in  tiie  model.  Tlie  orcanization  of  this  report  will  follow  the  task  structure 
contained  in  the  original  proposal. 

During  the  course  of  the  work  carried  out  in  cooperation  with  the  AFGL, 
there  were  several  concurrent  developments  which  served  to  significantly 
advance  the  mesoscale  modeling  program.  The  first  of  these  events  was  the 
completion  of  a  thesis  by  D.  Medal  (1985),  which  takes  standard  surface  and 
uppe-'  air  observations,  interpolates  them  to  a  uniform  grid,  and  then  uses 
those  data  to  initialize  a  three-dimensional,  time-dependent  meso-alpha  model 
with  a  UO  kilometer  mesh  size  covering  an  area  1600  kilometers  on  a  side. 

This  model  is  then  run  for  a  period  of  3  hours  to  provide  initial  conditions 
for  a  meso-beta  model  with  a  10  kilometer  grid  size  covering  an  area  250 
kilometers  on  a  side.  The  thesis  thus  contains  a  real  data  initialization 
procedure  and  a  comparison  of  model-predicted  rainfall  with  observed  data  from 
a  severe  precipitation  event  in  the  south  of  France. 

The  main  contribution  of  Medal's  thesis  was  the  development  of  a  complete 
set  of  procedures  to  inialize  the  model  of  Nickerson  and  Richard  (1981)  with 
real  data.  A  second  and  concurrent  development  consisted  of  the 
implementation  of  a  new  time  integration  procedure,  the  incorporation  of  an 
upper  wave  absorbing  layer,  and  the  testing  of  improved  lateral  boundary 
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conditions.  There  is  now  a  new  version  of  the  model  (Nickerson,  et  al . ,  1986) 
which  will  replace  the  version  of  the  meso-heta  model  used  by  Medal. 

The  nine-track  magnetic  tapes  whose  contents  are  described  in  Tables  2-4 
have  already  been  delivered  to  AFGL,  together  with  a  copy  of  Medal's  thesis, 
an  English  translation  of  a  part  of  that  thesis,  and  a  listing  (in  Fortran)  of 
all  the  programs  in  Medal's  code.  The  model  was  installed  at  the  computing 
center  at  Los  Alamos  National  Laboratory  and  required  little  in  the  way  of 
modification  to  the  existing  code.  All  of  the  subroutines  were  run  and  the 
output  then  compared  with  the  calculations  carried  out  in  France.  Installing 
the  model  code  at  AFGL  should  only  require  a  change  in  the  I/O  procedures. 


2.  TASK  I:  Real  Data  Initialization 

A  systematic  investigation  of  a  heavy  rainfall  event  has  been  carried  out 
in  order  to  compare  model  predictions  with  observations.  A  thesis  oy  Medal 
(1935)  contains  a  detailed  description  of  a  procedure  which  uses  standard 
radiosonde  and  surface  observations  to  generate  a  balanced  set  of  initial 
conditions  on  a  regular  grid  suitable  for  initializing  the  meso-alpha  ve'^sion 
of  the  meteorological  prediction  model.  The  meso-alpha  model  is  similar  to 
the  meso-beta  model,  except  that  the  equations  are  modified  to  include  a  map 
factor,  and  the  precipitation  processes  have  been  eliminated.  The  model  has 
been  integrated  over  a  4l  x  4’  horizontal  grid  using  a  grid  length  of  40 
kilometers.  There  are  15  vertical  levels  in  both  the  meso-alpha  and  mesc-teta 
versions  of  the  model.  An  English  translation  of  the  introductory  section  is 
included  with  Medal's  thesis  submitted  as  part  of  this  final  report. 

I  The  area  under  investigation  is  a  mountainous  region  in  the  Cevennes,  a 

part  of  France's  Massif  Central  where  an  extensive  rain-gauge  network  has  been 

'y 
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installed  in  order  to  provide  warnings  of  flash  flooding.  This  region  has 
historically  been  subjected  to  heavy  rainfall  events  with  attendant  loss  of 
life  and  extensive  property  damage.  Hourly  rainfall  data  are  available  to 
provide  not  only  an  evaluation  of  the  total  rainfall  predicted  by  the  model, 
but  also  to  provide  some  insight  on  the  ability  of  the  model  to  depict  the 
temporal  evolution  of  the  storm. 

The  following  two  paragraphs  contain  a  translation  of  the  conclusions 
from  Medal's  thesis: 


The  rainfall  episodes  in  the  Cevennes  are  characterized 
by  a  transgression  of  hot  and  humid  Mediterranean  air.  The 
synoptic  configuration  and  the  orographic  effects  seem  to  be 
the  main  ingredients  giving  rise  to  rainfall;  this  was  the 
case  on  29  August  1976,  a  period  we  have  chosen  to 
simulate.  For  this  purpose,  we  have  developed  and  tested 
different  initialization  procedures  for  the  Nickerson  and 
Richard  model  (1981).  The  procedure  of  dynamic  adaptation 
of  the  meteorological  fields  to  the  sub-synoptic  scale, 
followed  by  a  spatial  nesting  of  the  mesoscale  model, 
produced  the  most  encouraging  results:  the  rainfall  field 
calculated  in  this  manner  displays  a  spatial  distribution 
which  corresponds  to  the  observations.  As  far  as  the 
simulated  quantities  are  concerned,  they  appear  to  be 
underestimated . 

Following  various  tests,  the  rainfall  calculation 
appeared  to  be  particularly  sensitive  to  the  initialization 
of  the  mass  field,  especially  temperature  and  humidity. 
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Perky  (1979)  had  already  stressed  the  influence  of  these  two 
fields  on  the  amount  of  simulated  rainfall.  Anthes  and 
Haagensen  (1983)  have  explained,  during  the  simulation  of 
catastrophic  rainfall  in  Sichuan  (July  12-15,  1979)  the 
deficit  of  calculated  versus  observed  rainfall  by  an 
underevaluation  of  the  initial  humidity  field  (Nickerson  et 
al.,  1984).  Our  own  conclusions  confirm  these  results  on  a 
smaller  scale  (10  km).  The  temperature  and  humidity  fields 
sensitive  to  local  factors  have  characteristic  scales  which 
are  probably  not  very  well  represented  by  the  mesh  of  the 
available  recording  network.  In  order  to  increase  the 
resolution  of  these  initial  fields  in  the  surface  layer,  a 
possible  solution  consists  of  introducing  the  climatological 
network  data  (pressure,  temperature,  humidity)  in  the 
initialization  procedure  on  numerical  models.  This  is  the 
subject  of  Appendix  III  which  remains  to  be  applied. 

Medal's  code  has  been  installed  and  run  on  a  Cray  1  at  Los  Alamos 
National  Laboratory.  Table  1  contains  a  simplified  flow  chart  showing  the 
organization  of  the  twelve  separate  programs.  A  listing  of  each  program, 
together  with  a  translation  of  the  comment  lines  is  included  with  this 
report.  The  codes  are  written  in  modular  form,  and  should  be  easily  adapted 
to  any  machine  having  a  Fortran  compiler.  Other  than  minor  modifications 
required  by  a  different  operating  system,  it  was  not  difficult  to  make  the 
codes  operational  on  the  Los  Alamos  computing  system.  Magnetic  tape  numbers 
9245  and  9223  contain  the  original  codes  and  sample  output  from  Medal's 
thesis.  The  third  tape  contains  the  codes  modified  to  run  at  Los  Alamos.  The 
contents  of  all  three  tapes  are  given  in  Tables  2-4. 
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During  the  period  6-9  November  1982,  the  Cevennes  was  once  again  visited 
by  an  exceptionally  severe  storm,  with  heavy  rains  accompanied  by  surface 
winds  in  excess  of  50  m/sec.  Captain  D.  V.  Ridge  of  AFGL  provided  us  with  3-D 
Neph  data  for  the  period  enabling  us  to  construct  cloud  cross-sections  through 
the  region.  Figure  1.1  shows  the  location  of  the  3-D  Neph  gridded  data  (the 
larger  grid),  and  also  the  inset  grid  showing  the  domain  covered  by  the  meso- 
beta  model.  The  line  labeled  NW  -  SE  indicates  the  location  of  the  schematic 
vo’^tical  cloud  cross  section  shown  in  Figure  1.2.  Figure  1.3  shows  the 
0'"irinal  working  drawing,  together  with  the  hand-written  data  garnered  from 
the  3“D  Neph  data.  Of  particular  interest  is  the  vertical  extent  of  the  cloud 
and  the  presence  of  rainfall  upwind  of  the  mountain  barrier  indicated  by  a 
surface  present  weather  code  of  "6." 

Originally  it  was  our  intention  to  directly  utilize  the  satellite  data  to 
initialize  the  meso-beta  moisture  fields  in  the  manner  described  in  an  earlier 
quarterly  progress  report.  However,  experience  gained  from  Medal's  thesis  led 
us  to  conclude  that  without  a  proper  balance  and  initialization  of  the  mass 
and  momentum  fields,  the  addition  of  satellite  data  would  do  little  to  improve 
the  forecast.  We,  therefore,  decided  to  use  a  radiosonde  located  in  the 
southeast  of  France  ('^3.5‘’N,  5‘’E)  together  with  cloud  data  constructed  from  12 
satellite  de'^ived  cross-sections  similar  to  Figure  1.2  to  initialize  the  meso- 
alpha  model. 

The  meso-alpha  model  was  initialized  with  a  single  sounding  corresponding 
t*j  ‘200  GMT  on  1‘  I.'ovember  ’982  and  run  f or  a  period  of  3  hours.  Medal's 
interpolation  program  was  then  used  to  construct  an  initial  data  field  for  the 
meso-beta  model,  which  was  subsequently  run  for  a  period  of  6  hours.  The 


following  section  describes  the  results  of  the  numerical  simulation. 


2.1  Model  simulation  of  7  November  1982 

The  initial  soundings  of  temperature,  vapor  mixing  ratio,  cloud  water, 
and  wind  are  given  in  Table  5.  In  order  to  determine  the  structure  of  the 
cloud  system  over  the  meso-beta  domain,  12  vertical  cross  sections  were 
constructed  using  the  3~D  Neph  data.  These  showed  that  cloud  bases  were  at 
approximately  500  meters,  and  that  cloud  tops  were  general 'y  5000  meters, 
except  over  the  higher  mountain  elevations  where  tops  were  at  approximately 
12000  meters. 

Two  runs  were  made  with  the  meso-alpha  model  (dx  »  90  km).  In  the  first 
simulation  the  cloud  water  mixing  ratio  was  set  to  2.0  g/kgm  in  the  height 
range  from  500  to  5,000  meters.  However,  when  the  surface  elevation  was  more 
than  1000  meters,  the  initial  cloud  extended  to  a  height  of  12,000  meters. 

For  heights  between  5,000  and  12,000  meters,  the  cloud  water  mixing  ratio  was 
set  to  1.0  g/kgm.  Whenever  cloud  water  was  initially  present  in  the  model, 
the  vapor  mixing  ratio  was  set  to  the  saturated  value  at  that  temperature.  A 
second  control  run  was  then  made  in  which  everywhere  in  the  model  the  cloud 
water  mixing  ratio  was  set  to  zero. 

A  comparison  of  the  two  simulations  showed  that  after  a  period  of  3  hours 
cf  model  integration,  the  two  simulations  were  essentially  identical.  That 
is,  the  presence  of  cloud  water  at  the  start  of  the  simulation  had  little 
influence  on  the  final  predicted  values. 

The  model  terrain  used  in  the  meso-Deta  calculation  is  shown  in  Figure  2, 
and  is  the  same  found  in  the  thesis  by  Vedal.  Surface  elevations  vary  from 
sea  level  to  a  height  of  1,700  meters.  The  Rhone  valley  separates  the 
Tevennes  in  the  northwest  from  the  Alps  in  the  east. 

Figures  3,1  through  3-7  show  the  evolution  of  the  wind  field  at  level  13 
during  the  course  of  the  6-hour  simulation.  One  must  remember  that  the  model 


makes  use  of  a  terrain-following  coordinate  system  and  that  this,  therefore, 
represents  the  flow  on  a  surface  which  varies  with  altitude.  Also,  the 
terrain  used  in  the  meso-beta  simulation  with  10  kilometer  grid  resolution  is 
not  the  same  as  that  used  in  the  meso-alpha  simulation  with  a  **0  kilometer 
grid  resolution.  It  takes  several  hours,  therefore,  for  the  airflow  to  reach 
a  quasi-steady  state  adjustment  to  the  new  underlying  terrain.  Horizontal 
plots  of  wind  speed  at  the  same  level  (Figures  ^.1  through  U.7)  show  that  the 
average  speed  c''ange3  very  slowly  at  the  end  of  the  6-hour  period  of 
integration.  The  adjustment  of  the  horizontal  wind  field  to  a  quasi-steady 
state  is  also  evident  in  the  plots  of  vertical  velocity  at  level  13  shown  in 
figures  6.1  through  5.6. 

The  initial  oloud  field  entirely  covers  the  model  domain.  However, 
Figures  6,1  through  6.7  show  that  the  higher  clouds  are  located  over  the 
Cevennes  and  the  Alps,  in  remarkable  agreement  with  the  conditions  reported  in 
the  3"D  Neph  data.  Nevertheless,  during  the  course  of  the  meso-beta 
simulation,  the  vertical  extent  of  the  clouds  decreases  substantially. 

Figures  7.1  through  7.7  show  the  cloud  water  mixing  ratio  at  level  13. 
There  is  a  very  rapid  adjustment  in  the  cloud  field  giving  way  to  maximum 
values  of  cloud  water  of  approximately  0.3  g/kgm.  It  should  be  recalled  that 
the  meso-alpha  model  does  not  contain  any  mechanism  for  a  conversion  of  cloud 
water  to  rain  water,  and  that  by  the  end  of  the  simulation  it  produces  large 
amounts  of  cloud  water.  During  the  first  hour  of  the  meso-beta  simulation, 
■"air.  drops  formed  and  grew  very  quickly,  thereby  depleting  the  initial  cloud 
water.  The  corresponding  plots  of  rain-water  mixing  ratio  are  shown  in 
Figures  8.1  through  8.6,  The  average  grid  point  value  decreases  from  .8^1  to 
,29  g/kgm  during  the  first  hour  and  then  decreases  at  a  much  slower  rate.  One 
especially  encouraging  aspect  of  this  simulation  is  that  clouds  and 


precipitation  are  produced  over  the  relatively  flat  terrain  upwind  of  the 
major  terrain  forcing.  Previous  runs  with  the  meso-beta  model  using  a  single 
sounding  initialization  procedure  yielded  only  orographically  forced  clouds 
over  the  elevated  terrain. 

In  some  locations,  observed  precipitation  rates  amounted  to  nearly  10 
mm/hr  for  the  3”day  storm.  Maximum  precipitation  rates  predicted  by  the  model 
were  in  excess  of  7.5  mm/hr,  although  a  detailed  analysis  of  the  observed 
rainfall  is  not  available  for  inclusion  in  this  report.  Tne  3-D  Neph  data 
indicated  a  present  weather  code  of  "6"  over  most  of  the  model  domain,  thereby 
indicating  the  presence  of  rain  and  lending  support  for  the  general  rainy 
conditions  in  the  numerical  simulation  (Figures  9.1  through  9.6). 

One  of  the  major  problems  experienced  by  AFGL  users  of  the  meso-beta 
model  has  been  the  appearance  of  significant  noise  in  the  pressure  field. 
Figures  10.1  through  10.3  show  the  surface  pressure  during  the  last  3  hours  of 
the  simulation.  It  is  seen  that  they  essentially  mirror  the  terrain  shown  in 
Figure  2,  and  do  not  give  any  indication  of  pronounced  oscillations  at  the 
boundaries . 

In  order  to  examine  the  three-dimensional  structure  of  the  simulation,  we 
now  show  several  vertical  cross  sections.  Figures  11.1  through  11.7  show  the 
"U"  component  of  the  horizontal  wind,  and  Figures  12.1  through  12.7  the  "V" 
component  along  a  west-east  cross  section  at  a  distance  of  170  kilometers  from 
the  southern  boundary.  These  figures  confirm  that  the  winds  at  upper  levels 
also  reach  a  quasi-steady  state  by  the  end  of  the  period  of  integration.  The 
same  also  holds  for  the  vertical  velocity  shown  in  Figures  13.1  through  13.6, 
although  the  zero-contour  line  shows  more  variability. 

The  previously  mentioned  large  initial  values  of  cloud  water  coupled  with 
zero  rain  water  lead  to  a  rapid  change  in  the  cloud  water  mixing  ratio  as 


shown  in  Figures  14, i  through  14.7.  It  takes  approximately  1  hour  for  the 
clouds  to  reach  a  state  which  subsequently  changes  by  rather  small  amounts. 

It  should  be  kept  in  mind,  however,  that  the  smallest  contour  is  0.1  g/kgm, 
and  that  the  actual  cloud  boundaries  extend  beyond  that  shown  by  the  0.1 

contour , 

The  evolution  of  the  rain  water  mixing  ratio  shown  in  Figures  15. i 
through  15.6  provides  further  evidence  of  the  microphysical  chain  of  events. 
Whereas  the  rain  water  was  initially  zero,  at  the  end  of  1  hour  the  maximum 
value  is  in  excess  of  3.0  g/kgm.  Here  again  the  minimum  contour  value  is  0.1 
g/kgm,  and  the  extent  of  the  rain  "cloud"  is  actually  larger  than  that  shown 
by  the  0.1  contour.  The  persistence  of  the  tall  cloud  along  the  eastern 
boundary  where  the  model  domain  extends  to  the  Alps  is  probably  due  in  part  to 
the  constancy  of  the  thermodynamic  variables  on  the  inflow  boundary.  The 
relatively  small  values  of  rain-drop  concentration  obtained  during  the  initial 
phase  of  the  simulation  are  due  to  the  unrealistic  initialization  of  the  cloud 
and  rain  parameters  (Figures  16.1  through  16.6).  The  initial  cloud  field 
■'esulted  in  an  almost  explosive  growth  of  large  rain  drops  during  the  first 
few  minutes,  followed  by  a  change  to  a  more  normal  growth  regime. 

Figures  17  through  19  show  selected  vertical  cross  sections  after  3  hours 
of  model  time  for  distances  of  100  km,  150  km,  and  200km  from  the  southern 
boundary.  Figure  17.4  is  especially  noteworthy  in  that  it  shows  the  presence 
of  ''ain  at  the  low  elevations  of  the  Rhone  Valley  at  this  cross  section  lOO  km 
from  the  southern  boundary. 

"Tnere  will  likely  be  a  follow-up  simulation  of  this  storm  by  the  modeling 
group  in  Clermont-Ferrand,  France,  where  they  are  attempting  to  improve  the 
initialization  procedure  for  mass  and  momentum  in  order  to  more  accurately 
relate  the  satellite-observed  moisture  fields  and  the  initial  model 
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/  parameters.  The  determination  of  vertical-velocity  fields  which  are 

S 

p  consistent  with  the  cloud  fields,  seems  especially  important.  Any  forthcoming 

.'  reports  on  this  subject  by  the  Laboratoire  Associ6  de  M§t6orologie  Physique 

(LAMP)  in  Clermont-Ferrand  will  be  forwarded  to  AFGL. 
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Figure  1.1 


Tlic  grid  .ire. I  represents  the  region  of  3“D  Neph  data.  The  inset  box  Is 
the  model  donuiine  for  the  meso-beta  Simulation,  and  the  line  'NW-SE" 
slu'ws  the  loc.it  ion  of  cross  sections  in  Figures  1.2  and  1.3. 
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Figure  1.3 

OriRiii.il  1-D  Noph  dnt.i  used  to  construct  Figure  1.2. 
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Figure  3.1 

Initial  horizontal  wind  field  for  sigma  level  13 
obtained  from  meso-alpha  simulation.  Minimum  and 
maximum  values  of  wind  speed  are  also  given. 
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Figure  3.2 

Same  as  Figure  3.1  but  after  1  hour. 
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Figure  3.3 

Same  as  Figure  3.1  but  after  2  hours. 
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Figure  3.4 

Same  as  Figure  3.1  but  after  3  hours. 
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Figure  3.5 

Same  as  Figure  3.1  but  after  4  hours. 
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Figure  3.6 

Same  as  Figure  3.1  but  after  5  hours. 
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Figure  3.7 

Same  as  Figure  3.1  but  after  6  hours. 
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Figure  4.5 

Same  as  Figure  4.1  but  after  4  hours. 
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Figure  5.2 

Vertical  velocity  for  level  13  after  2  hours. 
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Figure  5.3 

Vertical  velocity  for  level  13  after  3  hours. 
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Figure  5.4 

Vertical  velocity  for  level  13  after  4  hours. 
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Figure  5.5 

Vertical  velocity  for  level  13  after  5  hours. 
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5.6 

Vertical  velocity  for  level  13  after  6  hours. 
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Figure  6.1 

Initial  vertically  Integrated  cloud  water 
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Figure  6.2 

Vertically  Integrated  cloud  water  after  1  hour. 


V- CLOUD 


(2:0:0) 


W-E  (KM) 


MIN  -+.000E+00  MRX  -+.133E+01 


V-CLOUD 


(3:0:0) 


H-E  (KM) 


MIN  -+.000E+O0  MRX  -+.130E-»-01 


MAX=  i  3049  XMAX=  250  YMAX=  90 

MiN=  0  XMIN=  160  YMIN=  130 

AVERAGE  VALUE=  285112130170 


Figure  6.4 

Vertically  Integrated  cloud  water  after  3  hours. 
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Figure  6.5 

Vertically  Integrated  cloud  water  after  4  hours. 
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Figure  6.6 

Vertically  Integrated  cloud  water  after  5  hours. 
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Figure  7.4 

Cloud  water  mixing  ratio  for  level  13 
after  3  hours. 
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Figure  7.6 

Cloud  water  mixing  ratio  for  level  13 
after  5  hours. 
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Figure  8.1 

Rain  water  mixing  ratio  for  level  13  after  1  hour. 
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Figure  8.3 

Rain  water  mixing  ratio  for  level  13  after  3  hours. 
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Figure  9.2 

Forecast  rainfall  after  2  hours. 
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Figure  9.5 

Forecast  rainfall  after  5  hours. 
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Figure  9.6 

Forecast  rainfall  after  6  hours. 
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Figure  10.1 

Surface  Pressure  after  4  hours. 
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Figure  10.3 

Surface  Pressure  after  6  hours. 
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Figure  12.1 

West-East  vertical  cross  section  of  V,  170  km  from  southern 
boundary  at  time  0. 
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Figure  12,3 

Same  as  Figure  12.1  but  after  2  hours. 
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Figure  12.4 

Same  as  Figure  12.1  but  after  3  hours. 
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Figure  12.5 

Same  as  Figure  12.1  but  after  4  hours. 
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Figure  13.3 

Same  as  Figure  13.1  but  after  3  hours. 
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Figure  13.6 

Same  as  Figure  13.1  but  after  6  hours. 
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Figure  14.6 

Same  as  Figure  14.1  but  after  5  hours. 
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Figure  15.1 

West-East  vertical  cross  section  of  rain  water  mixing  ratio, 
170  km  from  southern  boundary  after  1  hour. 
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Figure  15.3 

Same  as  Figure  15.1  but  after  3  hours. 
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Figure  15.4 

Same  as  Figure  15.1  but  after  4  hours. 
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Figure  16.3 

Same  as  Figure  16.1  but  after  3  hours. 
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Figure  16.5 

Same  as  Figure  16.1  but  after  5  hours. 


=  +.000E-»-00  MRX  =* 


4- .  6  2  7  E  4-  001 


Figure  16.6 

Same  as  Figure  16.1  but  after  6  hours. 
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Figure  17.3 

Same  as  Figure  17.1  but  for  cloud  water 
mixing  ratio. 
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Table  2 

Contents  of  Tape  Number  9245 


Format:  1600  bpi,  ASCII,  NL 

RECFM-FB 

RECORD  LENGTH- 80BL0CK  SIZE=400 


File  Number 


Program  Name 
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CRAY01-MEZM0DLS 

2 

CRAY01-MEZM0DMS 
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DOM.PRENU9 

5 
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6 

DOM.PRESOL2 
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DOM.PREBAL2 
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D0M.RW.INTBAR01 

9 
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10 
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1 1 
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Table  3 


Contents  of  Tape  Number  9223 


Format:  Same  as  tape  number  9245 

File  Number 


Program  Name 


1  CRAy.MZ3DLS.OUT05 
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3  CRAY.AMA3DRR 
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3.  TASK  II:  Time  Dependent  Lateral  Boundary  Conditions 

The  problem  of  specifying  lateral  boundary  conditions  in  the  mesoscale 
model  has  been  examined  in  some  detail,  as  part  of  the  general  upgrading  of 
the  numerical  procedures  carried  out  under  TASK  III  (see  section  ^  of  this 
report).  Unfortunately ,  however,  the  question  of  lateral  boundary  conditions, 
when  viewed  within  the  context  of  four-dimensional  data  assimilation,  has  not 
been  solved.  The  enclosed  conference  paper  by  Marroquin  and  Brown  (1985) 
clea'’ly  describes  a  fundamental  difficulty  in  properly  specifying 
mathematically  correct  lateral  bounda-^y  conditions  which  are  not  changing  in 
time.  They  conclude  their  paper  with  the  following  remarks. 

•'The  gravity  wave  experiments  reveal  that  there  exists  an 
incompatibility  between  the  absorbing  layer  diffusion  (second- 
order  or  Rayleigh)  and  the  lateral  radiation  boundary 
formulation.  This  incompatibility  is  manifested  in  reflection  of 
horizontally  propagating  gravity  waves  impinging  upon  the  lateral 
bo'inda’^ies  in  the  absorbing  layer  region.  Other  authors  have 
found  the  same  difficulty  (Van  der  Berg  and  Oerlemans,  1985).  We 
have  found  that  wave  reflection  is  substantially  reduced  by  using 
a  single  value  of  the  phase  velocity  for  the  boundary 
computation,  "his  value  is  the  average  of  the  phase  velocities 
computed  from  the  Orlanski  formulation  below  the  absorbing 
layer.  Modelers  should  be  aware  of  this  problem  that  seems  to  be 
more  prominent  for  short  gravity  waves." 

nevertheless,  the  real-data  initialization  procedure  reported  on  under  TASK  I 
above,  might  be  used  on  an  ad-hoc  basis:  The  meso-alpha  code,  MEZDLS,  (see 
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Table  1)  could  be  run  either  before  or  concurrently  with  the  meso-beta  code, 
MEZDMS.  At  every  hour  the  results  could  be  saved  as  input  for  the 
interpolation  code,  AMFMEZS.  This  new  set  of  initial  conditions  generated 
from  the  meso-alpha  model  would  then  be  incorporated  into  the  meso-beta  model 
using  a  weighting  function  which  has  the  value  of  unity  on  the  boundary,  and 
which  decreases  to  zero,  four  or  five  points  into  the  model  domain.  The  value 
on  the  boundary  would  then  contain  the  large-scale  initial  conditions  and  the 
value  some  five  points  into  the  meso-beta  model  would  contain  the  values 
predicted  by  that  model.  Those  points  in  between  would  serve  as  a  transition 


The  need  for  additional  theoretical  studies  of  this  problem  cannot  be 
overemphasized.  Moreover,  the  presence  of  underlying  terrain  further 


complicates  the  situation.  Vertical  coordinate  transformations  such  as  the 
sigma-P  system  used  in  the  present  version  of  the  model  offer  certain 


advantages,  but  they  also  add  an  extra  degree  of  complexity  to  the  problem. 


TASK  III:  Integration  Procedures 

A  new  version  of  the  model  has  been  developed  in  which  the  Shapiro  filter 
has  been  eliminated  and  replaced  with  an  explicit  horizontal  diffusion 
operator,  and  the  TASU  Matsuno  time-integration  procedure  replaced  by  an 
Asselin  filter.  The  new  version  of  the  model  is  described  in  the  accompanying 
paper  by  Nickerson,  et  al.,  (1986)  to  be  published  by  the  Monthly  Weather 
Review,  February  1986.  In  order  to  Implement  the  new  scheme  and  verify  tnat 
it  worked  properly,  it  was  found  necessary  to  conduct  an  extensive  series  of 
tests  using  a  two-dimensional  version  of  the  model. 


The  new  model  code  is  contained  in  the  third  file  on  tape  number  9223 
previously  given  to  AFGL.  For  integration  periods  of  3  to  6  hours  over  the 
Alsace  terrain  used  in  that  three-dimensional  simulation,  the  incompatibility 
problem  between  the  upper  wave  abso'^bing  layer  and  the  lateral  boundary 
condition  did  not  seem  to  be  especia’ly  serious.  Nevertheless,  such 
conclusions  might  not  hold  for  other  regions.  Additional  tests  should  be 
conducted  using  different  initial  conditions  over  the  same  terrain.  However, 
what  would  contribute  most  to  the  development  of  a  validated  model  is  the 
specification  of  a  set  of  well-posed  three-dimensional  numerical  experiments, 
the  '■esults  of  which  could  be  checked  against  independent  predictions  by 
either  theory  or  laboratory  experiment. 

5.  TASK  IV;  Aerosol  Model 

Considerable  progress  has  been  made  in  the  evaluation  of  the  aerosol 
model.  Although  three-dimensional  calculations  using  the  aerosol  model  have 
beer  made  by  Chaumerliac  et  al . ,  (1983),  an  extensive  series  of  tests  have 
just  been  completed  using  the  new  discretization  schemes  described  above  under 
TASK  in.  The  new  mic^ophysical  scheme  of  Richard,  et  al . ,  (1989)  now  permits 
the  aerosol  model  to  include  nucleation  scavenging,  a  significant  advancement 
over  the  previous  version  of  the  model. 

An  extensive  series  of  tests  have  been  carried  out  seeking  to  quantify 
the  role  of  the  various  physical  processes  in  the  aerosol  model.  For  example, 
the  scavenging  efficiencies  of  maritime  and  continental  clouds  are 
different.  An  article  is  presently  in  preparation  for  submission  to  the 
Jcurral  of  Geophysical  Research . 


Although  the  aerosol  model  has  not  been  extended  to  include  the  larger 
size  category  appropriate  for  the  treatment  of  dust  particles,  in  principle 
that  could  be  accomplished  by  adding  another  log-normal  distribution  similar 
to  one  presently  used  for  the  accumulation  mode.  In  that  regard,  J.  F. 

Kafouf,  a  student  at  the  LAMP  has  completed  a  thesis  containing  a  new 
formulation  for  the  lower  boundary  conditions.  The  new  PBL  formulation 
includes  a  separate  equation  for  the  turbulent  kinetic  energy,  and  the 
treatment  of  vegetation  and  soil  moisture.  The  ability  to  simulate  the  proper 
generating  conditions  for  dust  at  the  lower  boundary  then  provides  the 
oppo>'tunity  for  a  realistic  assessment  of  the  resuspension  and  deposition  of 
dust  in  the  atmosphere.  A  copy  of  that  thesis  will  be  sent  to  AFGL  when  it 
becomes  available. 

6.  CONCLUSIONS 

A  simulation  of  a  heavy  rainfall  event  has  been  carried  out  over  southern 
France  using  3-D  Neph  data  supplied  by  AFGL  to  initialize  a  meso-alpha 
model.  The  model  of  Medal  (1985)  was  run  for  a  period  of  3  hours,  and  the 
output  from  that  simulation  was  subsequently  interpolated  to  a  meso-beta  grid 
in  order  to  provide  balanced  initial  conditions  for  a  regional  scale  model. 

An  examination  of  the  results  shows  that  the  meso-beta  model  is  capable  of 
simulating  many  features  of  an  orographically  forced  precipitation  event, 
including  temporal  and  spatial  information  on  the  amount  of  liquid  water 
present  in  the  atmosphere,  as  well  as  the  number  concentration  (and  hence  size 
distribution)  of  rain  drops. 

The  model  is  capable  of  being  relocated  to  other  geographical  areas  for 
the  purpose  of  simulating  clouds,  precipitation  and  airflow  over  complex 


orography.  Despite  the  presence  of  rather  severe  topographical  conditions  on 
the  eastern  boundary  of  the  meso-beta  model  which  included  part  of  the  Alps, 
no  serious  computational  dificulties  were  encountered  in  either  the  meso-alpha 
or  the  meso-beta  simulations. 

7.  SUGGESTIONS  FOR  FURTHER  WORK 

In  order  to  interface  with  other  data  sets  and  output  from  larger  scale 
models  such  as  the  Air  Force  global  model,  it  may  be  desirable  to  use  a 
projection  ether  than  the  Lambert  projection  used  by  Vedal  to  calculate  the 
map  factor  terms  in  the  model  equations.  In  addition,  the  initialization 
procedure  should  be  modified  to  make  it  easier  to  move  the  model  domain  to  a 
different  geographical  location. 

At  the  present  time,  the  model  is  not  capable  of  utilizing  all  of  the 
detailed  moisture  information  contained  in  the  3-D  Neph  data  sets. 

Simulations  have  shown  that  the  mere  inclusion  of  realistic  initial  moisture 
dist'^ibutions  will  have  little  effect  on  the  resultant  simulations  until  a 
technique  is  developed  to  derive  balanced  motion  and  divergence  fields  which 
are  consistent  and  compatible  with  the  condensed  water  present  in  the  air. 

Just  as  the  meso-alpha  model  was  used  to  provide  balanced  initial 
distributions  of  wind  and  moisture  for  the  meso-beta  model,  a  way  must  be 
developed  to  provide  balanced  initial  conditions  for  the  meso-alpha  model 
usin,q  the  output  from  a  larger  scale  model  as  well  as  data  from  remote  and  in 
situ  observing  platforms. 
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